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§1. Introduction 



Abelian monopoles play a key role in the dual superconductor mechanism of 
confinement ^ in non-abelian gauge theories. Abelian monopoles appear after the 
^ ■ so called abelian projection 2 \ According to the dual superconductor mechanism a 

£f) , condensation of abelian monopoles should give rise to the formation of an electric 

flux tube between the test quark and antiquark. Due to a non-zero string tension 
the quark and the antiquark are confined by a linear potential. This mechanism has 
been confirmed by many numerical simulations of the lattice gluodynamics 3 **' 4 ) which 
show that the abelian monopoles in the Maximal Abelian projection are responsible 
for the confinement. The SU(2) string tension is well described by the contribution 
of the abelian monopole currents 5 ); these currents satisfy the London equation for 
Q_i| a superconductor 6 )' 7 ) . In Fig. 1, taken from Ref. 7 \ the abelian monopole currents 

near the center of the flux tube formed by the quark-anti-quark pair are shown. It is 
seen that the monopoles wind around the center of the flux tube just as the Cooper 
pairs wind around the center of the Abrikosov string. In Fig. 2 taken from Ref. 8 ) we 
^ ■ show the dependence of the value of the monopole condensate <?* n ^ on f3 is shown. 

It is clearly seen that vanishes at the phase transition and it plays the role 
of the order parameters 9 ) . In Ref. 10 ) the effective lagrangian for monopoles was 
reconstructed from numerical data for monopole currents for 5(7(2) gluodynamics 
in the Maximal Abelian gauge. It occurs that this lagrangian corresponds to the 
Abelian Higgs model, the monopole are condensed in the classical string tension of 
the Abrikosov-Nielsen-Olesen string describes well the quantum string tension of the 
SU(2) gluodynamics. It means that the description of the gluodynamics at large 
distances in terms of the monopole variables can be very useful. 
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Fig. 1. The monopole currents around the string tube which is formed between the static quark 
and anti-quark in gluodynamics, Ref. 7 \ 



Below we give a short review of the recently obtained results; for the elementary 
introduction to the subject see Ref. 4 \ 
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Fig. 2. The dependence of on /3, the size of the lattice is L A ■ 4, L — > oo. 
Since the abelian monopoles play a dominant role in the confinement phenomena 

it is important to understand how do the abelian monopoles arise in the non-abelian 

gauge theories. This question is well understood in the lattice gauge theory while it 

deserves some additional study in the continuum theory. In Section 2 we discuss the 

origin of the abelian monopoles in the compact electrodynamics. 

In Sections 3 and 4 we discuss the procedure of abelian projection in the contin- 
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uum. As an example we consider the Polyakov Abelian gauge (Section 3) and the 
Maximal Abelian gauge (Section 4). 

In Section 5 we show that the abelian monopoles carry also an abelian electric 
charge, this fact means that the abelian monopoles are dyons. The dyonic nature of 
the abelian monopoles is a consequence of the complicated topological structure of 
non-abelian vacuum which contains self-dual non-abelian configurations. Due to the 
presence of such configurations the abelian monopoles become abelian dyons. 

It is possible that the abelian monopoles correspond to some non-abelian objects 
which populate the vacuum of gluodynamics. If there exists such correspondence 
one should observe a non-zero excess of energy density near the monopole world 
trajectories. Indeed, this excess had been found in numerical simulations of the 
lattice SU(2) gluodynamics n ) . In Section 6 we discuss this question in details. 

In Section 7 we discuss a new three-dimensional Maximal Abelian projection 
which is better than the usual Maximal Abelian projection from the physical point 
of view. 



§2. Siglular Gauge Transformations and Monopoles 

At first we discuss the question how the abelian monopoles appear in the £7(1) 
gauge theory. We use the lattice regularization, the action is a periodic function of 
the field strength tensor: 

Stance = E f(Re[U x ,^])= E S p i aq .(F x ^); U x ^ = e lF ^ , (2-1) 

where U x ^ v = U X ^U X+ ^ )V U^ +U ^U^ V and sum is over all lattice plaquettes. 

The obvious requirement is that for the weak fields the lattice action (2-1) is 
reduced to the continuum action: 

Splaq.\p xtiv _Q— > F x ,^v (2"2) 

Since the function e lx is periodic, the field-strength F x ^ v is physically equivalent to 
Fx^u + 1im Xjllv , where n x ,\xv is an arbitrary integer-valued lattice two-form. 

In the continuum limit n x ^ u is a surface (5-function. Thus in the continuum 
theory the gauge fields which correspond to 

F^x) and to F^(x) + 2Tre^ a p E ( 2 ' 3 ) 

i 

are physically equivalent. Here Zjfy(x,x^) is an arbitrary surface J-function*** : 

4^(x,x«) = jd?af u 5^(x - x«(a«)), 

dVg = e^daxfd^tfo-W K ' 

the coordinates of the surface, x$(a^), are parametrized by aa\a = 1, 2. 
*^ This surface is nothing but the Dirac sheet. 
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Thus denned £7(1) gauge theory is invariant under a large class of the singular 
gauge transformations, — > A ll — if2 + d fl f2, where Q = e ia is a regular function (an 
element of the £7(1) group) while the function a may contain discontinuities a — > a + 
2-7T. Under singular gauge transformations, F M1 ,(.A) — ► ^^(^4 ) = F^ U (A) + [<? M , 
It may be shown that [d^, d v ]a corresponds to the surface 5-function, eq.(2-4). Hence 
the singular gauge transformation shifts the field strength tensor as follows: 

F^(x)^F^(x) + 2ne^ o 0^E®(x,x®). (2-5) 

i 

Thus S[F flu (A)} = S[F^(A Q )\, where S[F^(A)} is the action of the compact £7(1) 
gauge theory. The explicit example of such singular gauge transformation is the 
following: let a be the azimuthal angle in the polar coordinate system: 

xi±ix2 = rsmje ±ta , X3 = rcos7 (2-6) 

F^{x)[A Q ] = F^(x)[A] + [d^d u }a = 

(2-7) 

= F^(x)[A] + 27^,1^,2 - S lti2 S v ,i)6(xi)S(x 2 ) 

The singular part in the field-strength appears due to the topological non-triviality 
of the gauge matrix, but the action which corresponds to A^ and to A^ is the same. 

Now we discuss how the abelian monopoles appear in the path integral for- 
malism. In the continuum limit the compact electrodynamics is described by the 
following partition function: 

Z = jvAj Dzexpj-^ / d 4 x(V^(A) + 2vr^ a ^(x,x^ j , (2-8) 

this theory is manifestly invariant under the singular gauge transformations which 
shift the field strength tensor F^ u according to eq.(2-5). The integration in eq.(2-8) 
is over collection of all open and closed surfaces U. The simplest (but not unique) 
measure of the integration over the surfaces £ is the string integration measure 
described in Ref. 12 ) . 

It is simple to rewrite the partition function (2-8) in the monopole representation. 
First we introduce the additional antisymmetric tensor field G^ v and represent the 
partition function as follows: 

Z = JvaJvgJ Pxexpj- J d'x^jG 2 ^ 

+iG fM „(F f , u (A) +2vre^Z' a/ 3(x,x))^| . (2-9) 

Integrating over the field A^ we get a constraint: d^G^ = which can be resolved by 
the use of the regular field B^: G^ u = (2-K)^ 1 e IMUa /3d a Bp. Changing the integration 
in eq.(2-9) from G to B we get: 

Z = JvbJ Vjexp^-J a 4 x^F* u (B) +ij^x,x)B^y (2-10) 
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where g = 4ir/e is the charge of the monopole, and the vector x^ parametrizes the 
boundaries of the Dirac sheets E /Jtl/ (x,x): 



3^1 



d v E®{x,x). (2-11) 



The boundary of the Dirac sheet is the world trajectory of the monopole. The field 
in eq.(2T0) plays the role of the dual gauge field. 
Applying the Bardakci-Samuel formula 13 ) to eq.(2-10) we get: 

Z = JvbJ P<f>exp{-/ d 4 x^F^( J B) + i|(^ + ^)#| 2 )|,(2.12) 

where is the (complex) monopole field which carries the magnetic charge 1. 

The situation with the SU{2) (SU(N)) lattice gauge theory is similar. There 
are (singular in the continuum) gauge transformations which do not change the 
lattice action but create a string-like singularity in the continuum limit. We discuss 
the physical consequences of this fact in the separate publication, below we give an 
explicit example of such gauge transformation. 

Consider 



Q{x) 



cos \ — sin ia 
sin ^e ia cos | 



(2-13) 



where a and 7 are respectively azimuthal and polar angles in the polar coordinate 
system (2-6). The action of the lattice SU(2) gauge model is invariant under this 
gauge transformation. But in the continuum, 

f^a") = n+P^Aw - in+id^d^n = q+f^a)^ 

-2vr(*5 M ,i^ 2 - S tl , 2 5 Utl )a 3 5(x 1 )5(x 2 )0(-x 3 ) (2-14) 

Thus the continuum limit of the lattice action is different from the naive continuum 
limit, limSiattice / / d 4 x^Tr[F^ u ], due to the presence of the singular fields. 



§3. Abelian Projection 



For the sake of simplicity we explain the Abelian projection for SU(2) gauge 
group, the generalization to SU(N) gauge group is straightforward. 

The Abelian Projection suggested by 't Hooft 2 ) is a partial gauge fixing defined 
by the conditions of diagonalization of some functional X[A] with respect to gauge 
transformations f2, the functional X[A] belongs to the adjoint representation of the 
SU(2) gauge group: 

X[A] -» X[A^] = f2 + X[A]f2 . (3-1) 

In the abelian projection the matrix X [A] is diagonal and the theory possesses the 
U(l) gauge symmetry. The generator of this symmetry is the Cartan group generator 
o" 3 /2, where a 3 is the Pauli matrix. The diagonal nonabelian gauge field ^4 3 behaves 
as an abelian gauge field with respect to the residual U(l) gauge transformations 
(An * A^ iQ+^d^Quii), Q u{ i) = exp{icr 3 a}). 



6 M.N.Chernodub, F.V.Gubarev, M.I.Polikarpov and A.I.Veselov 



The abelian monopoles appear as additional abelian degrees of freedom which are 
associated with the singularities in the gauge fixing conditions. These singularities 
appear at the points of the space where eigenvalues of the matrix X[A] coincide. 
The proof that these singularities are monopoles is given in Ref. 2 ) . Two eigenvalues 
of the matrix X[A] coincide if three independent equations are satisfied 2 ) and in the 
four-dimensional space these singularities form closed loops*). 

Consider as an example the Polyakov Abelian projection in SU(2) gluodynamics. 
This abelian projection is defined for a finite temperature gauge theory or for the 
gauge theory in a finite box. The Polyakov Abelian projection corresponds to the 
diagonalization of the functional P x related to the Polyakov loop, P = ^TYP^: 



P x = Texp{i£ dx A (x)J , (3-2) 



where the integration is over the closed path C x which starts and ends in the same 
point x and is parallel to the "time" direction; the symbol T means the path ordering. 
The path is closed due to the periodic boundary conditions. 

Now we show that in the continuum all the abelian monopoles are static in the 
Polyakov abelian gauge**). Consider a monopole trajectory which passes through 
some point x, thus the eigenvalues of the Polyakov loop coincide with each other in 
the point x. This means that the matrix L x belongs to the center of SU(2) group, 
Zj2'- TrL-r = ±2. Consider another point y which lies on the same Polyakov loop 
(this means that yi = x«, i = 1, 2, 3). TrL y = TrL x = ±2 and the eigenvalues of the 
matrix L y coincide since this matrix belongs to the center of SU{2) group. Thus, if 
the abelian monopole passes through the point x = (to,x) it also passes through all 
points y with the same spatial coordinates: y = (t, x) for all t. Thus in the Polyakov 
Abelian projection all abelian monopole trajectories are static. 

Now we discuss the Faddeev-Popov gauge fixing procedure for the 't Hooft 
abelian projection. The conditions of diagonalization can be explicitly written as 
follows: 

C a [A}=0, C a [A] =Ti{X[A]a a ) , a = 1,2, (3-3) 

where o~ a are the Pauli matrices. The Faddeev-Popov determinant ZVfpLA] is defined 
by the path integral over the gauge orbits of SU(2) group: 



a=l,2 

Straightforward evaluation of this integral yields: 



A FP [A] I VQ [] 6 (c a L4 (fi) ]) • (3-4) 

J 1 o 



zVfpLA] = const. expf^J d 4 xln|AiL4(x)] - A 2 L4(x)]|j 

= const. exp|2 J d 4 xln |/m{Ai[A(x)]}|| , (3-5) 



The closeness of the monopole loops reflects the conservation of the magnetic charge. An 
explicit proof of this fact is given in Ref. 14 ^ . 

**) The same conclusion was independently obtained by F.Lenz, private communication. 
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where A a [A(a;)], a = 1, 2 are the eigenvalues of the matrix A[A]. The Faddeev-Popov 
determinant (3-5) is explicitly gauge invariant. 
Substituting unity (3-4) into the path integral 

Z = JvA exp{-S[A}} . (3-6) 

and integrating over Q we get the partition function in the abelian gauge (3-3): 

Z g . f . = fvA exp{-S[A}}A FP [A] U S (C a [A\) . (3-7) 

J a=l,2 

§4. Maximal Abelian Projection 

The most interesting results on abelian monopoles were obtained in the Maximal 
Abelian (MaA) gauge. This gauge is defined by the maximization of the functional 

maxi2[i«] , R[A] = ~jd i x [(A*) 2 + (A 2 ) 2 ] , (4-1) 

The condition of a local extremum of the functional R is 

{d ll ±igA%)Ai = Q. (4-2) 

This condition (as well as the functional i?[A]) is invariant under the £7(1) gauge 
transformations, A M — > A M + d^a. The meaning of the MaA gauge is simple: by 
gauge transformations we make the gauge field A M as diagonal as possible. 
The Maximal Abelian gauge on the lattice is defined by the condition 14 ) : 

m^R[Up], R[U l ]=Y,Tr[a z U+a z U l }, l = {x,fi}. (4-3) 
" i 

This gauge condition corresponds to an abelian gauge, since R is invariant under the 
U(l) gauge transformations. 

Now we discuss the Faddeev-Popov gauge fixing procedure for the MaA projec- 
tion in the continuum. We define the Faddeev-Popov unity: 

1 = App[A; X] ■ j Vf2 exp{Ai?[A^]} , A -» +oo , (4-4) 

where App is the Faddeev-Popov determinant. We substitute the unity (4-4) in the 
partition function (3-6), shift the fields by the regular transformation 1? + : A — > A Q+ 
and use the gauge invariance of the Haar measure, the action and the Faddeev-Popov 
determinant are invariant under the regular gauge transformations. Thus we get the 
product of the volume of the gauge orbit, / Vf2, and the partition function in the 
fixed gauge: 

Z MaA = fvA expj-^ f d 4 x F 2 „[A] + XR[A]j A FP [A; A] . (4-5) 
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In the non-degenerate case the FP determinant can be represented in the form: 
A FP [A; A] = Det L 2M[A^ aA ] exp{-Ai?[A^ MaA ]j + . . . , (4-6) 

where Q^ aA = Q^ aA {A) is the regular gauge transformation which corresponds to 
a global maximum of the functional ii[A fi ], the dots correspond to the terms which 
are suppressed in the limit A — > oo; and 



M? y [A] 



xyL 1 dco a (x)dco b (y) 



(4-7) 



uj=0 



Q(uj) = exp{iu; a T a }, T a = cr a /2 are the generators of the gauge group, a a are the 
Pauli matrices. In the limit A — > +oo the region of the integration over the fields A 
reduces to region where the gauge fixing functional R is maximal, and therefore the 
partition function (4-5) can be rewritten as follows 15 ^: 

Z M aA = J VA^ V {-S(A)}Det\ ^M[A}^ r FMR [A] , (4-8) 

where Ffmr[A] is a characteristic function of the Fundamental Modular Region 16 ) 
for the MaA projection 15 ^: Ffmr[A] = 1 if the field A belongs to the Fundamental 
Modular Region (the global maximum of the functional R[A\) and Ffmr[A] = 
otherwise. 

Usually in the abelian projection the U(l) gauge invariant quantities [O) are 
considered. Below we derive the explicit expression for the SU (2) invariant quantity 
O which corresponds to O. The expectation value for the quantity O in the MaA 
gauge (4-1) is: 

< O >MaA = ^— [ VA exp{-S(A) + XR[A}} A FP [A; A] 0(A) . (4-9) 

Shifting the fields U —> U Q and integrating over Q both in the nominator and in the 
denominator of expression (4-9) we get: 

<0> MaA =<0>, 0(A) = Jvf2eMXR[A o ]} . (4-10) 

O is the SU (2) invariant operator. In the limit A — > +oo we can use the saddle point 
method to calculate O: 

£ DetkM\A n \ 0(A Q ) 

& W = - n(a) ; . ^ 

£ Det2M[A n(k) ] 
k=l 

where fi^ are the ^-degenerate global maxima of the functional i?[A fi ] with respect 
to the regular gauge transformations Q: R[A Q( ~ 3) \ = R[A Q( ~ k) ], j, k = 1, . . . , N. In 
the case of non-degenerate global maximum (N = 1), we get O(A) = 0(A^ (1> ). 
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§5. Abelian Monopoles Carry Electric Charge 

Consider a (anti-) self-dual configuration of the SU(2) gauge field: 

F^{A) = ±-e^ a pF a p(A) = ±F M1 , , (5-1) 

where F tlu {A) = d^A^ + i[An,A u ]. In the MaA projection the commutator term 
Tr((7 3 [.A^, A u ]) of the field strength tensor is suppressed, since the MaA pro- 
jection is defined 14 -* by the maximization condition (4-1). Therefore, in the MaA 
projection eq.(5T) yields 17 ^: /^(A) = d^A^ — d v A\ m ±f^ u (A). Thus, the abelian 
monopole currents must be accompanied by the electric currents: = d u f^ u (A) ~ 
ztd l/ f tMl/ (A) = ±J™. Thus, in the MaA projection the abelian monopoles are dyons 
for (anti) self-dual SU{2) field configurations 17 \ Below we show that in the real (not 
cooled) vacuum of lattice gluodynamics the abelian monopole currents are correlated 
with the electric currents 18 -* . 

In order to study the relation of electric and magnetic currents, we have to cal- 
culate connected correlators of these currents. The simplest correlator <C JJ?J^ S>= 
< JJTJu > — < JJT >< Ju > is zero, since < J™J^ >= due to the opposite 
parities of the operators J m and J e , and < J™' 6 >= due to the Lorentz invariance. 
The simplest non-trivial (normalized) correlator is 

G = -^<J?{y)Wv)q{y)>, (5-2) 

where q(x) is the sign of the topological charge density at the point x and p m>e = 
J2i < \ J] n ' e \ > /(4V) are the densities of the magnetic and the electric charges, V is 
the lattice volume (total number of sites) . 



f \ 

* i 

3» : 
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Fig. 3. The dependence of the correlator G on {3 

We perform the numerical calculation of the correlator (5-2) in the SU{2) lat- 
tice gauge theory on the 8 4 lattice with periodic boundary conditions. We use 100 
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statistically independent gauge field configurations for each value of 0. 

The dependence of the correlator G on (3 is shown in Fig. 3. This correlator is 
positive for all values of (5. Therefore, the abelian monopoles in the MaA projection 
carry an electric charge. According to definition (5-2), the sign of the electric charge 
of the monopole coincides with the product of the magnetic charge and the topologi- 
cal charge. Thus, in the gluodynamic vacuum the abelian monopoles become abelian 
dyons due to a non-trivial topological structure of the vacuum gauge fields 18 ) . 

§6. Abelian Monopole Currents are Correlated with SU (2) Action 

Density 

Abelian monopoles appear as singularities in the gauge transformations 2 ) " 4 ) 
(see also Sections 2-4). On the other hand, the monopole currents reproduce the 
SU(2) string tension 5 ). Thus, monopoles are likely to be related to some physical 
objects. A physical object is something which carries action. Below we study the 
local correlations of the abelian monopoles with the density of the magnetic and the 
electric parts of the SU(2) action (the global correlation was found in Ref. 19 )). We 
show that the monopoles are physical objects but it does not mean that they have 
to propagate in the Minkowsky space; a chain of instantons can produce a similar 
effect: an enhancement of the action density along a line in Euclidean space. The 
simplest quantities which can show this correlation are the relative excess of the 
magnetic and the electric action densities rj M,E = (S^ ,E — S)/S in the region near 
the monopole current. Here S is the expectation value of the lattice plaquette action, 
Sp =< (1 — T;TrUp) >. The quantities S^' E are, respectively, the magnetic and the 
electric parts of the SU(2) action density, which are calculated on plaquettes closest 
to the monopole current. 
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Fig. 4. The relative excess of the magnetic (triangles, from Refs. and the electric (boxes) action 
density near the monopole current. 
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Fig. 5. The relative excess of the magnetic action density near the monopole current for l 3 (squares) 
and 2 3 (triangles) monopoles on 24 4 lattice. 



In the continuum notation, the quantities S^ ,E have the following form: 
S m=\< TrMx) iV(x)) 2 > , S E = \ < Tr(n M (x) F^(x)f > , (6-1) 

n^x) is the unit vector in the direction of the current: n^(x) = j^(x)/\j^(x)\, if 
jfj,(x) ^ 0, and n^(x) = if j^(x) = 0. It is easy to see that for a static monopole 
(jo 7^ 0; ji = 0,i = 1,2,3) (Sm) corresponds to the chromomagnetic action 
density (Bf) 2 (chromoelectric action density (Ef ) 2 ) near the monopole current. 

We calculate the quantities t] M and i] E on the lattice 24 4 with periodic boundary 
conditions. In Fig. 4 we show the quantities rj M,E vs. [3. The monopole currents are 
calculated in the MaA projection. In Fig. 4 the statistical errors are smaller than 
the size of the symbols. It is clearly seen that the abelian monopoles are correlated 
with the magnetic and the electric parts of the SU(2) action density. Note that the 
correlation of the monopole current with the magnetic action density is larger than 
the correlation of the monopole current with the electric action density. 

The similar results are obtained with extended monopoles 20 ) which are defined 
on the cubes of the size N x N x N 21 \ In Fig. 5 we show the quantities i] M vs. (3 
for the extended monopoles of the sizes n 3 , n = 1,2 on the lattice 24 4 . 

§7. 3D Maximal Abelian Gauge and Effective Monopole Potential 

As we already discussed the largest part of the numerical calculations is per- 
formed in the MaA projection. Usually the expectation values of abelian operators 
(operators constructed from the abelian gauge fields, diagonal gluons) are calculated 
in this projection. But abelian operators correspond to nonlocal in time operators 
in terms of the original SU(N) fields U x ^ (see eq. (4-11)). This nonlocality occurs 
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Fig. 6. The effective potential for the monopole fields $ (for <P > 0) in the confinement (a) and the 
deconfinement (b) phases for the 3D MaA gauge. 



since the gauge fixing condition (4-1) contains time like links U x ^. For time nonlo- 
cal operators there are obvious problems with the transition from the Euclidean to 
Minkowsky space-time. Thus there are problems with physical interpretation of the 
results obtained for abelian operators in the MaA projection. 

Below we present first results of calculations in the 3D MaA projection which is 
defined by the same maximization condition as the usual MaA projection (4-1), but 
the summation in R (4-1) is over the space-like links*). Since the time-like links are 
excluded from the gauge fixing condition, the abelian operators in the 3D MaA pro- 
jection correspond to local in time operators constructed from the nonabelian fields. 
In ref. 8 ) the effective potential for the monopole creation operator was calculated in 
the MaA projection. As we discussed in the Introduction in the confinement region 
(below the critical temperature) this potential is of the Higgs type, above the criti- 
cal temperature this potential has minimum at the zero value of the monopole field. 
Also this behavior of the potential is very reasonable it is important to prove the 
monopole condensation in the 3D MaA gauge. In Fig. 6 (a,b) we show the effective 
monopole potential for the confinement phase ((3 = 1.5, 12 3 • 4 lattice) and for the 
deconfinement phase (/3 = 2.5, 12 3 • 4 lattice). It is clearly seen that the minimum of 
the potential is at the nonzero value of the monopole field for the confinement phase 
and is at the zero value for the deconfinement phase. Our definition of the effective 
potential V(<P) is the same as in Ref. 8 ): 



e -vm =<6 ^-^ mon (x))>, 



(7-1) 



here <& m(m {x) is the monopole creation operator, defined in ref. 8 ). We discuss the 
dependence of the minimum of the effective potential on the temperature for the 3D 
MaA projection in a separate publication. 

This gauge was discussed by U.-J. Wiese in 1990, was recently rediscovered by D. Zwanzinger 
(private communication to M.I. P.), and discussed by M. Miiller-Preussker at this school. 
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